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ABSTRACT 

Context. Theories on the origin of magnetic fields in massive stars remain poorly developed, because the properties 
of their magnetic field as function of stellar parameters could not yet be investigated. Additional observations are of 
I utmost importance to constrain the conditions that are conducive to magnetic fields and to determine first trends about 
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their occurrence rate and field strength distribution. 

Aims. To investigate whether magnetic fields in massive stars are ubiquitous or appear only in stars with a specific 
' spectral classification, certain ages, or in a special environment, we acquired 67 new spectropolarimetric observations 

^ ll for 30 massive stars. Among the observed sample, roughly one third of the stars are probable members of clusters at 

I , different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. 

O ■ Methods. Spectropolarimetric observations were obtained during four different nights using the low-resolution spec- 

tropolarimetric mode of FORS 2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of 
^ ' the VLT. Furthermore, we present a number of follow-up observations carried out with the high-resolution spectropo- 

^ ' larimeters SOFIN mounted at the Nordic Optical Telescope (NOT) and HARPS mounted at the ESO 3.6 m between 

2008 and 2011. To assess the membership in open clusters and associations, we used astrometric catalogues with the 
highest quality kinematic and photometric data currently available. 
' Results. The presence of a magnetic field is confirmed in nine stars previously observed with FORS 1/2: HD 36879, 

^ ! HD 47839, CPD-28 2561, CPD-47 2963, HD 93843, HD 148937, HD 149757, HD 328856, and HD 164794. New mag- 
. netic field detections at a significance level of at least 3a were achieved in five stars: HD 92206c, HD 93521, HD 93632, 

' CPD— 46 8221, and HD 157857. Among the stars with a detected magnetic field, five stars belong to open clusters with 

' high membership probability. According to previous kinematic studies, five magnetic O-type stars in our sample are 

• , candidate runaway stars. 

^ , 

, Key words, polarization - stars: early-type - stars: magnetic field - stars: massive - stars: kinematics and dynamics — 

. open clusters and associations: general 

1. Introduction the origin and the role of magnetic fields in massive stars 

>< During the last years a gradually increasing number of O, , 

H , 1 TTTTi X 1 ^ ■ i • ^ 1 r i he present study IS the contmuation oi our magnetic 

. early B-type, and WR stars have been mvestigated tor mag- -ii, t?/-voc i /o i- j-u ^rrrr c i r 

. . . „ .-'.^ ' 1, , , , ,. ^ , surveys with FORSl/2 at the VLT of a sample ol mas- 

netic fields, and as a result, about a dozen magnetic O-type - Ajoi i - li 4-j-a- i 

,^ riTii i ill sive O and B-type stars m open clusters at dmerent ages 

and two dozens ot early B-type stars are presently known j • ii, c i j , , j i J ■ nnnr; / tt • j- i 

/ T-i i 1 i ^AA^ i Tj 1 • J- 1 lonno oniT-n l orn lu i and m the field Started already m 2005 (e.g., Hubrig et al. 

(e.g., Donati et al. Hubrig et al. [2008, 201^ |2Qllbl r^oj^ j 2011 we presented the results of 41 snec- 

Martins et al.^OTU', Wade et al.[2nnilMl2a, 2012b). The re- i"™^ \ Vn^ a o P''™!^'^ tne results ot 4t spec 

^ , ^ ^. „ ^. „ — — . ' ^ tropolarimetric l(jRb2 observations ot 36 massive stars, 

cent detections ot magnetic helds m massive stars generate . ^i i i i ii i ir r^i ^ 

^. ^. ^'^ ^ , ^, , ^. , ^ Among the observed sample, roughly hall ot the stars were 

a strong motivation to study the correlations between evo- , , r i ^ ^ j-a- ^ i, • -t i 

, ° . , ., in .,. probable members oi clusters at ditterent ages (Hubrig et al. 

lutionary state, rotation velocity, and surtace composition ir,m i u h rpi • t, j r ■ r iT 

,„ , rrjTrrn t-. j. i i ^ah^ k i ^ j ^ j I2011bl) . 1 he new survey IS based On tour nights ot observa- 

(Brott et al. 120111 Potter et al. 120121) . and to understand 7: -.i r^r^Tioo • onn ix^ ^i, i. 

^ ' tions with 10Rh2 m 2011 May. I'urthermore, we present 

* Based on observations obtained at the European Southern ^ number of follow-up observations carried out with the 

Observatory, Paranal, Chile (ESO programmes 087.D-0049(A) high-resolution spectropolarimeters SOFIN mounted at the 

and 187.D-0917(A)), and SOFIN observations at the 2.56m Nordic Optical Telescope (NOT) and HARPS mounted at 

Nordic Optical Telescope on La Palma, the ESO 3.6 m between 2008 and 2011. The information 

** F.R.S.-FNRS Postdoctoral Researcher, Belgium on the occurrence rate of magnetic fields and field strength 
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Table 1. List of massive stars observed for this pro- 
gramme. Spectral classifications are listed according to the 
Galactic O Star Catalogue (Mai'z-Apellaniz et al. 2004|. 



Name 


Other 


V 


Spectral Type 




Identifier 






HD 36879 


BD+21899 


7.58 


05 HI (f) var 


HD 47839 


15Mon 


4.64 


06.5 HI (n)(f) 


CPD-28 2561 


CD-28 5104 


10.09 


06.5 f?p 


CPD-47 2963 


CD-474551 


8.45 


04 HI (f) 


HD 77581 


Vela X-1 


6.93 


B0.5 I ae* 


HD 92206c 


CD-57 3378 


9.07 


08.5 V p 


HD 303225 


CD-59 3225 


9.74 


B* 


HD 93026 


CD-58 3518 


9.67 


B2 HI* 


HD 93521 


BD+38 2179 


7.04 


09.5 V p 


HD 93632 


CD-59 3328 


8.38 


05 III (f) var 


CPD-58 2611 


CD-58 3526 


9.55 


06 V ((f)) 


HD 93843 


CD-55 5846 


7.33 


05 III (f) var 


HD 130298 


CD-55 5846 


9.29 


06.5 HI (n)(f) 


HD 148937 


CD-47 10855 


6.77 


06.5 f?p 


HD 149757 


C Oph 


2.58 


09.5 V nn 


HD 328856 


CD-46 11016 


8.50 


09.5 HI* 


CPD-46 8221 


CD-46 11017 


9.19 


+* 


HD 150958 


CD-46 11019 


7.29 


06.5 la (n)f+ 


HD 151804 


HR6245 


5.25 


08 laf 


HD 153426 


CD-38 11431 


7.49 


09 II-III 


HD 153919 


CD-37 11206 


6.53 


06.5 la f+ 


HD 154643 


CD-34 11503 


7.15 


09.5 m 


HD 156041 


CD-35 11432 


9.2 


BO V* 


HD 156154 


CD-35 11445 


8.04 


08 lab (f) 


HD 157857 


BD- 10 4493 


7.81 


06.5 HI (f) 


HD 164794 


9Sgr 


5.93 


04 V ((f)) 


HD 170580 


BD+03 3727 


6.69 


B2 V* 


HD 170783 


BD+04 3778 


7.73 


B5* 


HD 172175 


BD-07 4642 


9.44 


06 I (n)f 


HD 191612 


BD+35 3995 


7.84 


08 f?p 



Notes. All objects were observed with FORS 2, except 
HD 36879, HD 47839, and HD 191612, which were observed with 
SOFIN, HD 151804, which was observed with HARPS, and 
HD 164794, which was observed with HARPS and SOFIN. 
* indicates spectral types taken from SIMBAD. 



distribution is critical for answering the principal question 
of the possible origin of magnetic fields in massive stars. 

In the following, we present 67 new measurements of 
magnetic fields in 30 massive stars. Our observations and 
the obtained results are described in Sects. [5] and [31 and 
their discussion is presented in Sect. E) 

2. Observations and results 

2.1. FORS 2 observations 

The spectropolarimetric observations were carried out in 
2011 May 4-8 in visitor mode at the European Southern 
Observatory with FORS 2 mounted on the 8-m Antu 
telescope of the VLT. This multi-mode instrument is 
equipped with polarisation analyzing optics, comprising 
super-achromatic half-wave and quarter-wave phase re- 
tarder plates, and a Wollaston prism with a beam diver- 
gence of 22" in standard resolution modeQ- Polarimetric 
spectra were obtained with the GRISM 600B and the nar- 
rowest slit width of 0'.'4 to achieve a spectral resolving power 

^ The spectropolarimetric capabilities of FORS 1 were moved 
to FORS 2 in 2009. 



of i? ^ 2000. The use of the mosaic detector made of blue 
optimised E2V chips and a pixel size of 15 /im allowed us 
to cover a large spectral range, from 3250 to 6215 A, which 
includes all hydrogen Balmer lines from H/3 to the Balmer 
jump. The spectral types and the visual magnitudes of the 
studied stars are listed in Table [TJ 

A detailed description of the assessment of the longitu- 
dinal magnetic-field measurements using FORS 2 was pre- 
sented in our previous papers (e.g., Hubrig et al. I2004al 
I2004b[ and references therein) . The mean longitudinal mag- 
netic field, (Bz), was derived using 

V _ gefteA^ 1 d/ ^^^^ ^ 

where V is the Stokes parameter that measures the circular 
polarisation, / is the intensity in the unpolarised spectrum, 
QeB is the effective Lande factor, e is the electron charge, 
A is the wavelength, nie the electron mass, c the speed of 
light, d//dA is the derivative of Stokes /, and (B^) is the 
mean longitudinal magnetic field. Errors are determined in 
the original spectra from photon statistics and propagated 
to the resulting V/I spectra. The error of {B^) is obtained 
from the formal error of the error-weighted least-squares 
fit. 

Longitudinal magnetic fields were measured in two 
ways: using only the absorption hydrogen Balmer lines or 
the entire spectrum including all available absorption lines. 
The lines that show evidence of emission were not used to 
determine the magnetic field strength. In addition, a recti- 
fication was carried out for a few massive stars with contin- 
uum slopes detected in their Stokes V spectra. The feasibil- 
ity of longitudinal magnetic field measurements in massive 
stars using low-resolution spectropolarimetric observations 
has been demonstrated by previous studies of O and B-type 
stars (e.g., Hubrig et al. ,2006. ,200gl I2009al 1201 lal 1201 Ibp . 

Ten targets were observed once, 13 were observed two 
or three times. HD 148937 was observed four times to assess 
the magnetic field variability over the rotation cycle. Apart 
from this star, which has a rotation period of seven days 
(Naze et al. I2010|) . no exact rotation periods are known for 
the other stars in our sample. The star ( Oph was observed 
nine times with the aim to search for magnetic/rotational 
periodicity (see Sect. [3]). 

The results of our magnetic field measurements are pre- 
sented in Table [5J In the first two columns, we provide the 
star names and the modified Julian dates at the middle of 
the exposures, followed by the signal-to-noise ratio in the 
individual spectra in Col. 3. In Cols. 4 and 5, we present 
the longitudinal magnetic field (i?z)aii using the whole spec- 
trum, and the longitudinal magnetic field (-Bz)hyd using 
only the hydrogen lines. All quoted errors are la uncertain- 
ties. In Col. 6 we indicate the wavelength regions used for 
the magnetic field measurements. The spectral regions con- 
taining the interstellar lines of the Ca ii doublet are usually 
not used in the measurements. In Col. 7 we list the rotation 
phases for the Of ?p star HD 148937, which is the only star 
in our FORS 2 observations with a known period of about 
seven days. Finally, in the last column we identify new de- 
tections by ND, previous detections by PD and confirmed 
detections by CD. 

In our sample of 25 stars observed with FORS 2 in 2011 
May the presence of a magnetic field is confirmed in six stars 
previously observed with FORS 1/2 (Hubrig et al. [MM 
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Fig. 1. Normalised FORS 1/2 Stokes / spectra of 0-type stars with detected magnetic fields observed from 2005 to 2011. 
Well known spectral lines are indicated; all Balmer lines from the Balmer jump to H/3 are visible. The spectra were offset 
from 1 by multiples of 0.5 for clarity. 



I2011bl I2011c|) : CPD-282561, CPD-472963, HD93843, 
HD 148937, C Oph, and HD 328856. New magnetic field de- 
tections were achieved in five stars: HD 92206c, HD 93521, 
HD 93632, CPD-46 8221, and HD 157857. Although pre- 
vious, mostly single observations of the five 0-type 
stars HD 130298, HD 153426, HD 153919, HD 154643, and 
HD 156154 indicated the possible presence of weak mag- 
netic fields (Hubrig et al. ,2011b,) . our recent observations 
do not reveal magnetic fields at a significance level of 
3(T in these stars. We note that these recent observations 
have been carried out on single epochs for four out of five 
stars, only HD 130298 was observed on two different nights. 
Since the rotation periods for them are unknown, the non- 
detection of their magnetic fields can probably be explained 
by the strong dependence of the longitudinal magnetic field 
on the rotational aspect. 

The strongest magnetic fields are detected in the two 
Of?p stars CPD-28 2561 and HD 148937. Walborn ([W5)) 
introduced the Of?p category for massive O stars dis- 
playing recurrent spectral variations in certain spectral 
lines, sharp emission or P Cygni profiles in He i and the 
Balmer lines, and strong C ill emission lines around 4650 A. 
Only five Galactic Of?p stars are presently known: HD 108, 
NGC 1624-2, CPD-28 2561, HD 148937, and HD 191612 



(Walborn et al. I2010p . and all of them show evidence for 
the presence of magnetic fields (Martin et al. 120101 Wade 
et al. 2012a, Hubrig et al. 2008, 2011b, Donati et ah [2006)) . 
The record holder is the faintest star in the sample of Of ?p 
stars, NGC 1624-2, for which a longitudinal magnetic field 
of the order of 5 kC was recently detected by Wade et al. 

dMUD. 



Besides the small group of Of?p stars, our measure- 
ments indicate that the presence of magnetic fields can be 
expected in stars of different classification categories and 
at different evolutionary stages, from young main-sequence 
stars up to the supergiant stages. The different spectral 
appearance of stars with confirmed and new detections is 
illustrated in Fig. [1] where we present normalised FORS 1/2 
spectra in integral light. Among the whole sample, the two 
Of?p stars CPD-28 2561 and HD 148937 exhibit the most 
conspicuous spectral variability clearly visible in Fig. (TJ 
showing the spectra of these stars obtained at two differ- 
ent epochs. To calculate the rotation phase of the observed 
spectra of HD 148937, we used the rotation period of 7.03 d 
determined by Naze et al. (20Q8J- The stars with confirmed 
and new magnetic field detections are discussed in more 
detail in Sect. |31 
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Table 2. Longitudinal magnetic fields measured with F0RS2 in the studied sample. Entries related to previous mea- 
surements are indicated by an asterisk in the second column of the Table. The wavelength region used for the magnetic 
field measurements is indicated in the sixth column. The seventh column provides the rotation phases for the Of?p star 
HD 148937, which has a known rotation period of about seven days (Naze et al. I2010p . In the last column ND stands 
for new detections, PD for previous detections (Hubrig et al. 1200811201 Ibp . and CD for confirmed detections. All quoted 
errors are la uncertainties. 



Name 


MJD 


S/N 


[G] 


(-^2) hydr 

[G] 


Spectral 
Regions (A) 


Phase 


Comment 


CPD-28 2561 


*55338.969 


2262 


-381±122 


-534±167 


3580-4588 




PD 




55685.982 


2481 


99±82 


65±108 










55686.984 


2516 


-44±80 


6±102 










55687.980 


2490 


269±81 


281±90 






CD 


CPD-47 2963 


* 55337.094 


2722 


-190±62 


-154±96 


3696-5622 




PD 




55686.019 


2691 


-66±57 


-119±82 










55687.037 


3035 


-169±52 


-77±76 






CD 




55688.013 


2164 


95±69 


148±101 








VelaX-1 


55686.096 


3604 


-15±30 


-11±45 


3562-5886 








55687.058 


2894 


-80±32 


-114±50 










55688.048 


2433 


57±37 


88±69 








HD 92206c 


55686.135 


2393 


204±46 


140±88 


3658-5944 




ND 




55687.085 


3439 


-104±30 


-123±54 










55688.151 


2528 


10±43 


23±46 








HD 303225 


55688.220 


1902 


-59±54 


-81±90 


3447-5299 






HD 93026 


55687.247 


1647 


-19±75 


-52±100 


3470-4900 






HD 93521 


55686.042 


3891 


32±52 


45±73 


3689-5414 








55687.015 


3304 


-202±66 


-293±91 






ND 




55688.063 


3059 


-74±62 


-64±68 








HD 93632 


55687.217 


2927 


-152±50 


-120±62 


3318-4593 




ND 




55688.198 


2443 


-10±32 


5±68 








CPD-58 2611 


55687.151 


2631 


-58±44 


-80±59 


3580-4588 








55688.125 


2005 


-25±56 


-9±74 








HD 93843 


*55339.099 


3504 


-157±42 


-173±56 


3370-5545 




PD 




55686.115 


3230 


-34±38 


-12±52 










55687.185 


3385 


-132±39 


-93±56 






CD 




55688.180 


2647 


-63±47 


-130±61 








HD 130298 


*55339.159 


3037 


113±38 


193±62 


3573-5284 




PD 




55686.347 


2264 


-64±66 


-64±85 










55688.261 


2699 


59±55 


122±72 








HD 148937 


*54550.416 


2616 


-276±88 


-145±104 


3210-3599 


0.631 


PD 




*55337.285 


3364 


-204±71 


-225±103 




0.529 






*55339.206 


2253 


-290±85 


-389±129 




0.802 


PD 




55686.370 


2664 


-110±90 


-165±115 




0.172 


CD 




55687.265 


3083 


-87±74 


-137±95 




0.299 


CD 




55688.285 


2941 


-139±33 


-174±121 




0.444 


CD 




55689.361 


2522 


-241±91 


-229±118 




0.597 


CD 



2.2. SOFIN and HARPS observations 

Four 0-type stars, HD 36879, 15Mon, 9Sgr, and 
HD 191612, previously obs erved with FORS 1 and SOFIN 
(Hubrig et al. i2008. 2009a[ l^giU)) . were foUowed-up by nine 
additional SOFIN spectropolarimetric observations from 
2008 September and 2010 July. Spectral classifications of 
these stars are presented in Table \T\ Among these stars, 
only for the Of?p star HD 191612 was the presence of a 
magnetic field of the order of a few hundred Gauss previ- 
ously reported by Donati et al. (|2006p . We used the low- 
resolution camera {R w 30 000) of the echelle spectrograph 
SOFIN (Tuominen et al. 1999) mounted at the Cassegrain 
focus of the Nordic Optical Telescope (NOT). With the 2K 
Loral CCD detector we registered 40 echelle orders partially 
covering the range from 3500 to 10 000 A with a length of 
the spectral orders of about 140 A at 5500 A. The polarime- 
ter is located in front of the entrance slit of the spectro- 
graph and consists of a fixed calcite beam splitter aligned 
along the slit and a rotating super-achromatic quarter-wave 



plate. Two spectra polarised in opposite sense are recorded 
simultaneously for each echelle order providing sufficient 
separation by the cross-dispersion prism. Two to four sub- 
exposures with the quarter-wave plate angles separated by 
90° were used to derive circularly polarised spectra. A de- 
tailed description of the SOFIN spectropolarimeter and its 
polarimetric data reduction is given in Ilyin (|2012p . 

The spectra were reduced with the 4A software pack- 
age fllvin l^nnn]) . Bias subtraction, master flat- field correc- 
tion, scattered light subtraction, and weighted extraction 
of spectral orders comprise the standard steps of spectrum 
processing. A ThAr spectral lamp is used for wavelength 
calibration, taken before and after each target exposure to 
minimise temporal variations in the spectrograph. The typ- 
ical signal-to-noise ratio (S/N) for observations of 9Sgr is 
between 280 and 380, while for HD 36879 and HD 191612 
it is about 200. The highest S/N ratio of 440 was achieved 
for 15Mon. 
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Name 


MJD 


C /AT 

b/N 


[G] 


/ R \ 

[G] 


bpectral Fnase 
Regions 


Comment 


C Oph 


* 54609.343 


4083 


105±45 


54±54 


3645-4960 


PD 




55686.321 


3719 


-164±62 


-218±70 




CD 




55687.116 


3239 


-83±63 


-90±81 








55687.287 


2659 


-153±75 


-222±97 








55687.385 


2592 


-101±88 


-84±104 








55688.103 


1964 


18±113 


18±144 








55688.240 


2792 


95±68 


97±80 








55688.335 


2226 


-72±94 


-37±108 








55688.433 


2486 


-16±89 


14±110 








55689.349 


2193 


58±90 


52±117 






HD 328856 


*55336.370 


3239 


-173±53 


-155±65 


3589-5669 


PD 




*55339.223 


2975 


-149±48 


-75±72 




PD 




55686.417 


2309 


-43±48 


-96±70 








55687.351 


2612 


9±46 


21±64 








55688.301 


2279 


-128±50 


-250±77 




CD 


CPD-46 8221 


55686.393 


2389 


-132±42 


-65±63 


3569-5158 


ND 




55688.319 


2195 


-83±45 


-69±68 






HD 150958 


55686.431 


1957 


133±65 


122±118 


3695-4974 






55687.372 


2206 


-105±63 


-74±106 






HD 153426 


*55336.338 


3446 


-27±53 


-10±62 


3650-4874 






* 55339.246 


3136 


-171±55 


-275±70 




PD 




55687.434 


2318 


19±65 


62±82 






HD 153919 


*55337.341 


2727 


-213±68 


~119±95 


3680-4452 


PD 




55687.421 


2727 


26±82 


152±102 






HD 154643 


* 55339.340 


3382 


110±34 


121±52 


3376-4910 


PD 




55689.374 


2925 


73±42 


163±58 






TT"pv 1 c^^n/i 1 


ooooy.oyo 


1 c^m 
iijUi 


1 iztoy 


1 1 no 


oDoy— OoDo 




HD 156154 


•55337.358 


2753 


-118±38 


-167±54 


3617-4974 


PD 




55688.418 


2498 


19±54 


117±72 






HD 157857 


55687.401 


3035 


-110±46 


-226±61 


3383-4974 


ND 




55688.348 


2785 


-116±50 


-51±72 






HD 170580 


55688.403 


2579 


-55±31 


-77±47 


3343-6196 




HD 170783 


55688.389 


3016 


-19±34 


-2±43 


3326-5885 




HD 172175 


55688.370 


2210 


-73±94 


-79±124 


3670-4480 





In addition, the two 0-type stars, HD 151804 and 9Sgr, 
previously analyzed using FORS 1 data (Hubrig et al. I2008P 
were observed on 2011 May 25 and 26 with a S/N of about 
280 with the HARPSpol polarinieter (Snik et aLf^DTTjl feed- 
ing the HARPS spectrometer (Mayor et al. I2003P at the 
ESO 3.6 m telescope on La Silla within the framework of the 
program 187.D-0917(A). This HARPS spectropolarimetric 
material became recently publically available in the ESO 
archive. The HARPS archive spectra cover the wavelength 
range 3780-6913 A, with a small gap around 5300 A. Each 
HARPSpol observation is usually split into four to eight 
sub-exposures, obtained with four different orientations of 
the quarter- wave retarder plate relative to the beam splitter 
of the circular polarimeter. The reduction was performed 
using the HARPS data reduction software available at the 
ESO headquarters in Germany. The Stokes / and V param- 
eters were derived following the ratio method described by 
Donati et al. (|1997p . ensuring in particular that all spuri- 
ous signatures are removed at first order. Null polarisation 
spectra (labeled with n in Table 121) were calculated by com- 
bining the sub-exposures in such a way that the polarisation 
cancels out, allowing us to verify that no spurious signals 
are present in the data. 

HARPS and SOFIN spectra have been measured using 
the moment technique previously introduced by Mathys 
who discussed at length the diagnostic potential of high- 
resolution circularly polarised spectra using this technique 



in numerous papers (e.g., Mathys 119931 I1995al I1995bp . 
Wavelength shifts between right- and left-hand side circu- 
larly polarised spectra are interpreted in terms of a longi- 
tudinal magnetic field (i?z)- {Bz) is obtained by perform- 
ing a least-squares fit, forced through the origin, of the 
wavelength shifts — Al, as a function of 2'gAX2. In this 
process, the individual measurements are weighted by the 
respective l/a^{Xu — Al). The uncertainty of the derived 
value of the longitudinal magnetic field is estimated by its 
standard error a{{Bz)), derived from the least-squares anal- 
ysis (e.g., Mathys & Hubrig [T995p . The results of the longi- 
tudinal magnetic field measurements using mostly He i and 
He II lines are presented in Table |31 

The high-resolution spectropolarimetric observations 
confirm our previous detections of magnetic fields in the 
stars HD 36879 and 9 Sgr. The magnetic nature of the 
slowly rotating Of?p star HD 191612 (P = 538d ) was for 
the first time reported by Donati et al. (|2006p who ob- 
served it over four nights in 2005 June using the ESPaDOnS 
spectropolarimeter installed on the 3.6-m Canada-France- 
Hawaii Telescope and announced the detection of a mean 
longitudinal field (B^) = -220 ± 38 G. Hubrig et al. ((20101) 
reported on a single measurement of the magnetic field of 
HD 191612, (Bz) = 450 ± 153 G, using SOFIN observations 
acquired in 2008 September at rotational phase 0.43. Wade 
et al. (|201ip presented 13 new field measurements of this 
star demonstrating that the magnetic data can be mod- 
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Table 3. Magnetic field measurements of five 0-type stars 
using SOFIN and HARPS observations. Already published 
measurements are marked by asterisks. All quoted errors 
are la uncertainties. 



MJD/HJD 


Phase 


(Sz>,n [G] 


[G] 


Instr. 


HD 36879 


* 54345.389 




180±52 




F0RS2 


2454723.737 




-125±78 


53±75 


SOFIN 


2455195.615 




-276±89 


-64±92 


SOFIN 


15 Moil 


* 54609.968 




134±52 




F0RS2 


2455201.553 




-162±37 


-62±42 


SOFIN 


HD 151804 


* 53476.369 




-151±90 




F0RS2 


*53571.025 




68±65 




F0RS2 


*53596.096 




82±46 




F0RS2 


55708.356 




-104±38 


-31±36 


HARPS 


9Sgr 


*53520.357 




-114±66 




F0RS2 


* 53594. 119 




211±57 




FORS2 


*53595.096 




-165±75 




FORS2 


2455074.387 




187±56 


-33±60 


SOFIN 


2455075.384 




242±74 


72±76 


SOFIN 


2455078.396 




238±59 


66±66 


SOFIN 


2455081.382 




-265±65 


-42±62 


SOFIN 


55707.323 




210±42 


39±41 


HARPS 


HD 191612 


*2454721.347 


0.43 


357±58 


67±65 


SOFIN 


2455395.592 


0.68 


- 144±34 


-39±36 


SOFIN 


2455407.583 


0.68 


-154±63 


-67±70 


SOFIN 




Fig. 2. Longitudinal magnetic field variation of the Of?p 
star HD 191612 according to the 537.6 d period deter- 
mined by Howarth et al. (j2007|l . Red symbols correspond 
to ESPaDOnS observations, while blue symbols are our 
SOFIN measurements (Hubrig et al. I^OTUl and Table E]). 
Note that the measurement errors for both ESPaDOnS and 
SOFIN observations are of similar order. 



eled as a periodic, sinusoidal signal with a period of 538 d 
inferred from spectroscopy. The authors mention in their 
work that the SOFIN observation of the longitudinal field 
of HD 191612, (B.^) = 450 ± 153 G, does not agree very 
well with the observed field variation, which corresponds 
to (an essentially) consistently negative longitudinal mag- 
netic field. To check the consistency of our measurements 
we re-analyzed the previous SOFIN observations using an 
extended list of blend free spectral lines with known atomic 
parameters. The obtained results, (B^) = 357 ± 58 G, in- 
dicates, however, that contrary to the results of Wade et 
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Fig. 3. Variable line profiles in the SOFIN spectra of 
HD 36879 obtained at two different epochs. Left panel: 
Spectral variability of the Feiii line profile at A6048 and 
the blend Feiii-I-Heii at A6078. Right panel: Spectral vari- 
ability of line profiles presenting blends of the C ill and O ii 
lines at A4647 and A4660. 



al. (|201ip . the magnetic field of HD 191612 is indeed rather 
strongly positive in the phase 0.43. Clearly, more magnetic 
field observations of this star are needed to properly char- 
acterise the variation of the longitudinal magnetic field and 
the magnetic field geometry. In Fig.[5]we present our SOFIN 
observations together with those published by Wade et al. 

mm . 

The supergiant HD 151804 was previously observed 
three times in 2005 with FORS 1 (Hubrig et al. 5005]) . No 
field detection at a significance level of 3a was achieved 
indicating that the field should be rather weak if present 
at all. Our measurement (B^) = — 104±38G based on 
HARPSpol observations confirms these previous results. 



3. Results for individual targets with detected 
magnetic fields 

3.1. HD 36879 

A magnetic field at a 3.5a significance level in HD 36879 was 
for the first time detected in 2007 using FORS 1 observa- 
tions (Hubrig et al. I2008p . Our most recent high- resolution 
SOFIN observations in 2009 December confirm the pre- 
vious detection and reveal a change of the field polarity. 
Walborn (|2006|) reports that this star shows unexplained 
spectral peculiarities and/or variations, among them pecu- 
liar narrow, variable Siiv emission lines at A1394 and A1403 
(e.g. Walborn & Panek I1984p . Distinct spectral variability 
is also discovered in Feiii and the blended lines Ciii/Oii 
using our two high-resolution SOFIN spectra obtained on 
two different epochs. The variations of Feiii and Ciii/Oii 
line profiles are presented in Fig. |31 Interestingly, Hei lines 
appear slightly split. This appearance presented in Fig. H] 
could possibly indicate the presence of He spots on the stel- 
lar surface. 

3.2. 15Mon 



Hubrig et al. (|2011bp showed that this star is a probable 
member of the young, rich open cluster NGC 2264 located 
in the Monoceros OB 1 association. From the spatial dis- 
tribution of Ha emission stars. Sung et al. (2008) have 
concluded that the young open cluster NGC 2264 can be 
subdivided into three star forming regions (SFRs): two ac- 
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Fig. 5. Distinct variability observed in the spectra of 
CPD-28 2561 in the Hei A4026 (top left), Hell A4686 (top 
right), H/3 (bottom left), and Hei A5876 (bottom right) 
lines. 



Fig. 4. Line profiles of He i in the spectra of HD 36879 ob- 
tained at two different epochs. From top left to bottom right 
we show the lines Hei A4713, A4922, A5015, and A5876. 



tive SFRs - the region around the bright emission nebula 
in the southwest of the brightest star S Mon and the region 
to the north of the famous cone nebula - and surrounding 
these an elliptical-shaped halo region. The study of Teixeira 
et al. (2012) indicates that NGC 2264 has undergone more 
than one star-forming event, where the unembedded, in- 
substantial inner disk population and embedded thick inner 
disk population appear to have formed in separate episodes. 
Markova et al. (|2004p consider 15 Mon as a Galactic 0-type 
star with a mass of 32 Mq and T^q = 37 500 K. Speckle 
measurements of 15 Mon indicated the presence of a com- 
panion of spectral type 09.5 V n (e.g. Gies et al. I1993p . 
In very high S/N spectra becomes the companion notice- 
able by very broad underlying components to the strong 
Hei lines. According to Kaper et al. (|1996p and Walborn 
(|2006p shows this star distinct peculiarities in the spectra, 
which could be typical for stars possessing magnetic fields. 
The X-ray luminosity based on Chandra observations is 
log(Lx) = 31.82 and the ratio of X-ray to bolometric lumi- 
nosity log(Lx/-Z^boi) = —6.63 is slightly above the average 
(Ramirez et al. 120041 Naze 1200'^ . as may be expected from 
a magnetic star. The high-resolution X-ray spectrum of 15 
Mon is typical for a star with its spectral type and, appar- 
ently, not especially hard (Walborn et al. I2009p . 

The previously measured magnetic field using the full 
FORS 1 spectrum coverage was of positive polarity at a 
2.6cr significance level (Hubrig et al. I2009al) . A weak mag- 
netic field was detected in this star at a significance level of 
4.4cr using SOFIN observations. No contribution of the sec- 
ondary component was detected in our FORS 1 and SOFIN 
spectra (see also Fig.[T]). 



3.3. CPD-282561 

Our previous observations of the Of?p star CPD— 28 2561 
enabled us to detect a magnetic field at the 3.1cr level using 
the whole spectrum and at the 3.2a level using Balmer 
lines (Hubrig et al. I2011b|) . Our new spectropolarimet- 
ric observations of this star, recently reported by us in a 
short IBVS article (Hubrig et al. I2012ap indicate that the 
magnetic field is variable. Due to the relative faintness of 
CPD— 28 2561 with ^=10.1, it was only scarcely studied in 
the past. Levato et al. (|1988p acquired radial velocities of 
35 OB stars with carbon, nitrogen, and oxygen anomalies 
and found variability of a few emission lines with a proba- 
ble period of 17 days. Walborn et al. (|2010p mentioned that 
CPD— 28 2561 undergoes extreme spectral transformations 
very similar to those of HD 191612, on a timescale of weeks, 
inferred from the variable emission intensity of the C ill 
A A 4647-4650-4652 triplet. Our spectropolarimetric obser- 
vations on three consecutive nights in 2011 reveal strong 
variations in several hydrogen and helium line profiles. The 
variation of Hei A4026, Hei A5876, and of H/3 on such a 
short time scale were not reported in the literature before. 
A few examples of the detected spectral line profile varia- 
tions are presented in Fig. [5j The contribution of a variable 
emission component is well visible in three lines. The pro- 
file of the Hei A4026 line appears only weakly variable. 
These variations have an amplitude of about 1-2% of the 
continuum level and can likely be explained by rotation 
modulation. The line Hell A4686 exhibited strong emis- 
sion during our observations in 2010, which transforms to a 
much weaker P Cygni line profile with a weak blue shifted 
absorption in 2011. A contribution of the emission compo- 
nent is clearly seen in the profile of the Hei A5876 line in 
2010, but its presence can only be suspected in the variable 
line depth observed in 2011. A strong emission component 
is well visible in the line center of H/3 in 2010. The shape of 
the line profile of H/3 is similar to the behaviour of the Ha 
line profile in the spectra of 19 Cep presented in Fig. 13 of 
Kaper et al. (1997J and in Fig. 1 of Kholtygin et al. ([2003)) . 

CPD— 28 2561 is an intrinsically bright X-ray source as 
detected by the Rosat X-ray observatory. Its ratio of X- 
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ray to bolometric luminosity (logLx/iboi ~ ^6) makes 
it unusually X-ray bright compared to other 0-type stars 
(Chlebowski [T^M Oskinova 2005). The hardness ratio 
(HRl 0.04 ±0.4) indicates that its X-ray emission can be 
harder than is typical in single 0-type stars. This provides 
indirect support for the presence of a magnetic field on the 
star, strong enough to confine the stellar wind (Babel & 
Montmerle .1997,) . 

3.4. CPD-472963 

For the star CPD— 472963, we confirm our previous S.lcr 
detection by measuring all absorption lines in the FORS 2 
spectra. The magnetic field of this star is variable and 
shows a change of polarity over our four-night observing 
run. According to Walborn et al. (P0TIJ)1 . CPD-47 2963 be- 
longs to the Ofc category, which consists of normal spectra 
with C III AA 4647-4650-4652 emission lines of comparable 
intensity to those of the Of defining lines N ill AA 4634-4640- 
4642. The origin of the magnetic field in this star probably 
differs from those of other magnetic 0-type stars, because 
non-thermal radio emission, which is frequently observed 
in binary systems with colliding winds, was detected by 
Benaglia et al. f2001). On the other hand, the membership 
of CPD— 472963 to a binary or multiple system has not yet 
been investigated. 

3.5. HD 92206c 

The presence of a variable magnetic field up to a 4.4(t signif- 
icance level in HD 92206c is documented here for the first 
time. This star is a probable member of the open clus- 
ter NGC3324 located inside of a partial ring of nebulos- 
ity northwest of the Carina Nebula. The assessment of the 
membership using astrometric catalogues with the highest 
quality kinematic and photometric data currently available 
is described in Sect. HI According to the Galactic O Star 
Catalogue (Mai'z-Apellaniz et al. belongs HD 92206c 

to the triple system HD 92206abc. An analysis of the distri- 
bution of the ionised and neutral material associated with 
the Hii region Gum 31 by Cappa et al. (2008) suggested 
that the excitation sources of this region are the brightest 
stars in the open cluster NGC 3324 with the triple system 
HD 92206abc considered as the brightest cluster member. 
The detection of protostellar candidates almost coincident 
in position with the open cluster NGC 3324 indicates that 
stellar formation is ongoing in this region. 

Walborn (fT982)) reports that the spectrum of HD 92206c 
displays very strong, broad hydrogen lines, possibly simi- 
lar to those in the Orion Trapezium cluster, in particular 
9^ OriC, and the appearance of these lines is indicative of 
extreme youth. A search for a magnetic field was carried 
out in another member of this cluster, the bright AO-type 
supergiant HD 92207. According to Hubrig et al. (j2012bp . 
this star also possesses a magnetic field. 

3.6. HD 93521 

A longitudinal magnetic field of the order of 200-300 G is 
detected in this star at a significance level of 3.2(t. This high 
Galactic latitude 09.5 Vp star with vsini — 390 km s^^ 
(Rauw et al. I2008P is one of the fastest rotators in our sam- 
ple, along with C, Oph. It is located at an unusually high 



Galactic latitude of 62.15 and displays prominent line pro- 
file variability in the optical and UV domain (Rauw et al. 
2008). In the most recent work by Rauw et al. (2012a|, the 
authors report that He and N are found to be overabundant 
and suggest that despite some ambiguities on the runaway 
status of the star it is likely that the star is a Population I 
massive 0-type star that was ejected from the Galactic 
plane either through dynamical interactions or as the re- 
sult of a supernova event in a binary system. 

HD 93521 is one of the few 0-type stars known to ex- 
hibit stellar pulsations. In this star, two frequencies of about 
2-3 hours are detected and interpreted as due to high- 
degree non-radial pulsation modes (Howarth & Reid 119931 
Rauw et al. 2008). 

The X-ray observations of this star have been summa- 
rized by Rauw et al. (|2012ap . They concluded that the stel- 
lar X-ray spectrum is consistent with a normal late 0-type 
star although with subsolar metallicity. No trace of a com- 
pact companion was found in the X-ray data. The ratio of 
X-ray to bolometric luminosity is log(Lx/iboi) — —7.1... — 
7.0 and the stellar X-ray luminosity Lx = 5 x 10^^ erg s~^ is 
rather low. Such low X-ray luminosities are also observed in 
other magnetic stars. Oskinova et al. (|201ip showed, that 
relatively low X-ray luminosities can be explained in the 
framework of the magnetically confined wind shock model 
if the stellar wind strength is correctly derived from an op- 
tical and UV spectroscopic analysis. Significantly, Rauw et 
al. (|2012ap show that besides a softer plasma component 
with kTi = 0.3 keV, which is typical for normal 0-stars, a 
harder component with kT2 = 3 keV is present in the X- 
ray spectrum of HD 93521. Recently, Ignace et al. (|2012p 
showed that such a hard component may serve as an indi- 
cator for the presence of a stellar magnetic field. 



3. 7. HD 93632 and HD 93843 

A weak negative magnetic field in HD 93632, just at a sig- 
nificance level of 3(7 was detected in one out of two FORS 2 
observations. Walborn (11973| reported that both HD 93632 
and HD 93843 appear slightly variable in their spectral 
types. Interestingly, also HD 93843, similar to HD 93632, 
exhibits a weak magnetic field. It's presence is confirmed in 
our new observations. According to Walborn et al. (|2010p . 
the star HD 93843 belongs to the Ofc category. Prinja 
et al. (|1998p monitored the stellar wind of this star us- 
ing lUE time series. They identified systematic changes 
in the absorption troughs of the Siiv and Nv resonance 
lines with a repeatability of wind structures with a pe- 
riod of 7.1 days. The authors suggested at that time that 
a magnetic field might help to explain the cyclical wind 
perturbation. Both stars have been observed during the 
Chandra X-ray survey of the Carina Nebula (Townsley et al. 
IMD)- For HD 93632 is the X-ray luminosity logix = 31.7. 
The ratio of X-ray to bolometric luminosity is small, only 
log(Lx/iboi) = —7.7. A one-temperature fit to the X-ray 
spectrum shows a slightly higher than average temperature 
kT = 0.56 (Naze et al. 20TTp. For HD 93843 is the X-ray 
luminosity logLx = 32.1. The ratio of X-ray to bolomet- 
ric luminosity is standard, log(Lx/^boi) = —7.1. A one- 
temperature fit to the X-ray spectrum derives a spectral 
temperature kT = 0.3, which is lower than in HD 93632 
(Naze et al. l^OTT)) . 
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Fig. 6. Longitudinal magnetic field variation of the Of ?p 
star HD 148937 according to the 7.032 d period deter- 
mined by Naze et al. (|2010p . Red symbols correspond to 
the ESPaDOnS observations (Wade et al. I2012b}, while 
green symbols are our FORS 1 and FORS 2 measurements 
(Hubrig et al. [MM [20TTbl and present paper). Note that 
the measurement errors for both ESPaDOnS and FORS 1 /2 
observations are of similar order. One more measurement 
(not shown in this Figure) was obtained in 2010 on May 22, 
but ordinary and extraordinary beams in the FORS 2 setup 
were overlapping in this exposure due to problems with 
slitlets. 



3.8. HD 148937 

The first detection of a mean longitudinal magnetic field 
(B^) = -254 ± 81 G in the Of?p star HD 14893 7 usin g 
FORS 1 at the VLT was reported by Hubrig et al. (pn05| . 
An extensive multiwavelength study of HD 148937 was car- 
ried out by Naze et al. (2008), who detected small-scale 
variations of Hen A4686 and the Balmer lines with a pe- 
riod of seven days and an overabundance of nitrogen by 
a factor of four compared to the Sun. The periodicity of 
spectral variations in hydrogen and helium lines was re- 
confirmed using additional higher resolution spectroscopic 
material indicating the similarity to the other Of?p stars 
HD 108 and HD 191612 (Naze et al.^UWI. Our new and old 
observations of HD 148937 are presented in Fig. [B] together 
with the ESPaDOnS observations obtained at the CFHT 
by Wade et al. ((2012bp . This figure demonstrates the excel- 
lent potential of FORS 2 for the detection and investigation 
of magnetic fields in massive stars: While an exposure time 
of 21.5 h at the CFHT was necessary to obtain seven binned 
measurements, the exposure time for FORS 2 observations 
accounted only for two to four minutes and only 2.3 h were 
used for our observations at six different epochs including 
telescope presets and the usual overheads for readout time 
and retarder waveplate rotation. 

3.9. C Oph 

The star Ophiuchi is a well-known rapidly rotating run- 
away star with extremely interesting characteristics. This 
star exhibits a wonderful bow shock as evident from IR 
Spitzer IRAC observations and undergoes episodic mass 
loss seen as emission in Ha. It is possible that ^ Oph rotates 
with almost break-up velocity with wsini = 400kms~^ 
(Kambe et al. I1993p . C Oph is the second target of our 
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Fig. 7. Phase diagram for the best sinusoidal fit corre- 
sponding to the period of 0.8 days for the longitudinal mag- 
netic field measurements using hydrogen lines for ^ Oph. 
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Fig. 8. Phase diagram for the best sinusoidal fit cor- 
responding to the period of 1.3 days for the longitudinal 
magnetic field measurements using the whole spectrum for 

COph. 



sample with detected non-radial pulsations with periods 
of hours, first discovered from ground-based observations 
(e.g. Kambe et al. 119971) . Later, the star was monitored by 
the MOST satellite, which revealed at least a dozen signifi- 
cant pulsation frequencies between 1 and 10 d~^, typical of 
/? Cep-like pulsation modes (Walker et al. I2005P . 

The first detection of a mean longitudinal magnetic field 
in this star was announced by Hubrig et al. (|2011cp . Nine 
additional FORS 2 spectropolarimetric observations show- 
ing a change of polarity have been obtained over all four 
observing nights. The rotation period of this star is still un- 
certain: Oskinova et al. pOOip studied ASCA observations 
that covered just the expected rotation period of the star 
of about 0.8-1.0 d. A X-ray flux variability with ~20% am- 
plitude in the ASCA passband (0.5-10 keV) with a period 
of 0.77 d was reported. The X-ray spectra of COph are not 
especially hard as found from Chandra HETGS and Suzaku 
data (Waldron & Cassineni i2007l Walborn et ah [2009]) . The 
variation period is also not well-defined in our observations. 
The resulting periodograms for the magnetic field measure- 
ments using hydrogen lines shows a dominating peak cor- 
responding to a period of about 0.8 d, while the measure- 
ments using the whole spectrum indicate a rotation period 
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Fig. 9. Combined IR Spitzer IRAC (3.6 /^m blue, 4.5 /xm 
green, 8.0 /xm red) image of the bow shock around VelaX- 
1. Archival data have been used (AOR 17858304, PI: R. 
Iping). Image size is 5' x 5'. 



of about 1.3 d. In Fig. [7] we present the magnetic field varia- 
tion using hydrogen lines phased with the period of 0.8 days, 
while Fig. |S] refers to magnetic field variations determined 
from all lines phased with 1.3 days. Interestingly, just a few 
months ago PoUmann (|2012p presented equivalent width 
changes of Hei A6678 with a period of 0.643 d, which is 
roughly half of the magnetic period of 1.3 d determined us- 
ing the whole spectral region. For the rotation period of 
~1.3d the equivalent width of the Hei line would display 
a double wave. Such a behaviour of Hei lines is frequently 
found in He variable early-type Bp stars (e.g. Briquet et al. 
[2004). Clearly, more measurements are needed to determine 
the true periodicity of the magnetic field in this star. 

A bow-shock nebula is also present close to the high 
mass X-ray binary VelaX-1, consisting of a pulsar in an 
eccentric (e — 0.09, Porb = 8.96 d) orbit around the here 
studied BO. 5 supergiant star HD 77581. In Fig.lHlwe present 
an image of the bow-shock nebula based on archival Spitzer 
IRAC maps (AOR 17858304, PI R. Iping). The mass-loss 
rate (M « lO~^M0yr~^) for VelaX-1 is well constrained 
from its accretion powered X-ray luminosity (e.g., Fiirst et 
al. [20T0l Oskinova et al. [20T2| . Contrary to our study of 
C Oph, no magnetic field detection at a significance level 
of 3(7 has been achieved in VelaX-1 during our observing 
run. However, we should keep in mind that the spectral 
behaviour of HD 77581 is much more complex than that 
of C Oph due to the presence of bumps and wiggles in the 
line profiles, and an impact of tidal effects producing or- 
bital phase-dependent variations in the line profiles lead- 
ing to asymmetries such as extended blue or red wings 
(Koenigsberger et al. 120121 . Furthermore, Kreykenbohm et 
al. (jMQSI detected flaring activity and temporary quasi- 
periodic oscillations in INTEGRAL X-ray observations. 



3.10. HD 328856 

Our new observations confirm the presence of a magnetic 
field in this star. As we already showed in our previous 
work on the basis of photometric membership probabil- 
ity, HD 328856 is a member of the compact open cluster 
Hogg 22 in the Ara region at an age of 5Myr and a dis- 
tance of about 1300 pc (Hubrig et al. 12011b]) . On the other 
hand, its proper motions indicate that HD 328856 is not 
fully co-moving with the other cluster members, deviating 
from the cluster mean proper motion by ^2a (for more de- 
tails of the membership probabilities, see Kharchenko et al. 

mm . 

3.11. CP 0-^46 8221 

A search for the presence of a magnetic field in this star 
was carried out for the first time. A weak mean longitudi- 
nal field was detected during the first observing night at a 
significance level of 3.1cr. As we show below in Sect. |4l this 
star has a high kinematical membership probability in the 
open cluster Hogg 22. 

3.12. HD 157857 

Spectropolarimetric observations of this star have been con- 
ducted for the first time. Our measurements on the second 
observing night carried out on hydrogen lines resulted in 
a field detection at a significance level of 3.7(7. Schilbach 
& Roeser ([2005]) retraced the orbit of HD 157857 in the 
Galactic potential and concluded that it was ejected from 
the open cluster NGC 6611 about 3.8 Myr ago with a pecu- 
liar velocity of 113.5kms~^. 

3.13. HD 164794 

All four SOFIN and one HARPS high-resolution spectropo- 
larimetric observations confirm our previous detection of 
the longitudinal magnetic field in this star. The field ap- 
pears to be variable with negative and positive extreme 
values of -265 G and -1-242 G. HD 164794 belongs to the 
young cluster NGC 6530 and according to the Galactic O 
Star Catalogue (Mai'z-Apellaniz et al. 12004^ . has a classifi- 
cation 04 V ((f)), exhibiting weak N ill emission and strong 
He II A4686 absorption. It is a well-known non-thermal ra- 
dio emitter and, according to van Loo et al. p006|) . the 
most likely mechanism is synchrotron emission from collid- 
ing winds, implying that all O stars with non-thermal radio 
emission should be members of binary or multiple systems. 
Hints of a wind-wind interaction were indeed detected in 
the X-ray domain (Rauw et al. I2002|l . A long-term study 
of its binary nature and spectrum variability has been re- 
cently presented by Rauw et al. (I2012bp who derived an 
orbital solution and an orbital period of 8.6 yr. 

3.14. HD 191612 

Our SOFIN spectropolarimetric observations of this star 
have already been discussed in Sect. 12.21 HD 191612 is cur- 
rently one of the best studied targets in the optical and 
X-ray domains among the five currently known Galactic 
Of?p stars. Due to the above mentioned discrepancies in 
the magnetic field measurements around the positive field 
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extremum, additional magnetic field observations are nec- 
essary for proper magnetic field geometry modeling. 

4. Discussion 

Whereas magnetic fields are ubiquitous in low mass stars 
due to their convective envelopes, and found in 10. ..20% of 
all intermediate mass stars, they were unexpected in O stars 
until just a few years ago. The first O star for which a mag- 
netic field was detected is the spectroscopically variable star 
9^ Ori C (Donati et al. 120021) . However, with ever more mag- 
netic O stars being identified in recent years (e.g., Hubrig 
et al. l2011b[ Martins et al. l2012p . it appears likely that the 
as yet small number of magnetic O stars is just a conse- 
quence of an observational bias, due to the low number of 
spectral lines in these objects. Indeed, there may be more 
reasons to find magnetic fields in O stars than in A and 
B stars. Fossil fields might appear with similar frequencies 
on both groups. If close binary effects as mass transfer or 
stellar merger would produce large scale stable surface fields 
in O stars (Ferrario et al. 2009] Langer [2012p . then about 
one quarter of them would be expected to be magnetic due 
to their very high close binary fraction (Sana et al. 120121) . 
Furthermore, their proximity to the Eddington limit pro- 
duces convection in the envelopes of O type stars, which 
could also produce surface magnetic fields (Cantiello et al. 

mm- 

With the present work we gradually increase the sample 
of massive stars with determined magnetic field strengths. 
It is the aim of our survey to provide new clues about the 
origin and detectability of the magnetic fields of massive 
stars. Bagnulo et al. (pUn|) used the ESO FORS 1 pipeline 
to reduce the full content of the FORS 1 archive. They claim 
that very small instrument fiexures, negligible in most of 
the instrument applications, may be responsible for some 
spurious magnetic field detections. We note however that 
out of the five Of?p stars, two were detected to be mag - 
netic for the first time by us (Hubrig et al. 12011bp . 
The excellent agreement of our FORS 1 and FORS 2 data 
with ESPaDOnS observations by Wade et al. (j2012b|) of 
HD 148937 can be seen in Fig. H Bagnulo et al. (POT^ 
could not detect the magnetic field in 9 Sgr (=HD 164794), 
while our high resolution observations of this star with 
both SOFIN and HARPS fully confirm our earlier result 
(Hubrig et al. l2008[) . Also, high resolution observations with 
ESPaDOnS (Petit et al. I2013[) showed a weak magnetic 
field in 40Tau (=HD 25558), which was already detected 
by Hubrig et al. P009bp . but not found by Bagnulo et al. 
(1201 2) from the same data. Overall, it is not clear why 
Bagnulo et al. (|2012p are not able to find magnetic fields 
found by other groups with FORS 1/2 and other instru- 
ments. These are only a few examples that show the need 
for confirmations of FORS 1 magnetic field measurements 
with high spectral resolution polarimeters. In the following, 
we discuss the role of magnetic stars in star clusters, mag- 
netic runaway stars, and indirect magnetic field indicators 
in massive stars. 

4.1. Magnetic stars in young clusters 

Of the 30 massive stars analysed in this paper, eleven stars 
have been found to be related to open clusters of differ- 
ent ages. The data on the cluster membership of these 
probable cluster 0-type stars are presented in Table El As 



database for the compilation of Table HI we used the All- 
sky Compiled Catalogue of 2.5 million stars (ASCC-2.5, 3rd 
version) of Kharchenko & Roeser (2009). The kinematic 
and photometric probabilities for cluster membership pre- 
sented in Cols. 3 and 4 were calculated according to the 
procedures described by Kharchenko et al. (|2004p . For cal- 
culating the kinematic membership probability, only proper 
motion information was used. Since the star HD 92206c 
in the open cluster NGC 3324 belongs to the triple sys- 
tem HD 92206abc, we present the proper motions of all 
triple system members in the same table. Largely deviat- 
ing proper motions of the cluster member HD 92206b in 
comparison to the average cluster proper motions are very 
likely caused by the orbital motion in the pair HD 92206ab. 
Interestingly, not only the cluster member HD 92206c shows 
the presence of a magnetic field. According to Hubrig et 
al. (2012b|), also another cluster member, the supergiant 
HD 92207, possesses a magnetic field of the order of a few 
hundred Gauss. Clearly, the cluster NGC 3324 should be 
included in future magnetic field surveys of massive stars 
to better understand the role of the environment on the 
generation of their magnetic fields. 

As we already mentioned in our previous work, ac- 
cording to Schilbach & Roeser ([20051) , the star HD 153919 
was probably ejected from the cluster NGC 6231 at an age 
of about 6.5 Myr. At present, HD 153919 is located 4.2° 
from the center of NGC 6231. Considering the proper mo- 
tions of HD 153919 given by Kharchenko & Roeser ([MM 
Ha = -1-1.42 ± 1.17mas/yr, fig = +4.84 ± 0.87 mas/yr) 
and those of the cluster by Kharchenko et al. (|2005a( 
-0.39 ± 0.52 mas/yr, -1.99 ± 0.38 mas/yr), then 2.15 Myr 
ago the star was at 0.55 ± 0.95° from the cluster center. 

Among the stars presented in Table only five stars 
show the presence of magnetic fields, HD 47839, HD 92206c, 
HD 328856, CPD-46 8221, and HD 164794. Two of them, 
HD 328856 and CPD-46 8221, belong to the open clus- 
ter Hogg 22, which appears as a second promising cluster 
for investigating the magnetic nature of cluster members. 
Although we deal here with small number statistics, we 
cannot refrain from mentioning that these five stars with 
detected magnetic fields belong to the youngest clusters in 
our small sample. 

4.2. Magnetic massive runaway stars 

Magnetic fields are detected in stars of spectral types, with 
different kinematical status, and at different evolutionary 
stages. Clearly, the detection of magnetic fields in all five 
Galactic Of?p stars implies a tight relation between the 
spectral characteristics of the Of ?p star group and the pres- 
ence of a magnetic field. Remarkably, apart from the Of?p 
star NGC 1624-2, no other Of?p star belongs to an open 
cluster or association. The analysis of the stellar content of 
the cluster NGC 1624 by Jose et al. ([^OTT]) indicated that 
NGC 1624-2 is the main ionizing source of the region, and 
the maximum post-main-sequence age of the cluster is es- 
timated as ~4Myr. 

Furthermore, recent studies of the evolutionary state of 
several magnetic massive stars indicate that some of them 
evolved significantly or have the status of a runaway star 
(e.g., Martins et al. [MB Hubrig et al. I2011b[ I2011dp . 
In the literature, the runaway status or candidate run- 
away status (i.e., if stars possess high space velocities 
determined through proper motion and/or radial veloc- 
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Table 4. Probable members in open clusters. 



Object 


ASCC 


Pkin 


-Pphot 


Cluster 


dist 


log* 




\ star 




My, <^t^y^ 


name 


number 


[%] 


[%] 




[pc] 


[yr] 


[mas/yr] 


[mas/yr] 


*HD 47839 


1021435 


62 


100 


NGC 2264 


660 


6.81 


-3.84, - 


'2.50,0.94 


-2.70, 


-3.50,0.25 


HD 92206a 


2231158 


80 


96 


NGC 3324 


2301 


6.72 


-10.32, 


6.76,4.05 


-7.81, 


3.77,0.64 


HD 92206b 


2231162 





12 


NGC 3324 


2301 


6.72 


49.61, 


0.86,3.70 


-7.81, 


3.77,0.64 


*HD 92206c 


2231147 


90 


97 


NGC 3324 


2301 


6.72 


-8.87, 


4.57, 1.90 


-7.81, 


3.77,0.64 


HD 93026 


2232306 


86 


88 


Bochum 10 


2027 


6.93 


-3.74, 


3.69,4.23 


-6.61, 


2.56,0.62 


CPD-58 2611 


2232406 


96 


100 


Trumpler 14 


2753 


6.67 


-3.09, 


4.44, 2.66 


-3.91, 


3.65,0.60 


*HD 328856 


2074567 


22 


100 


Hogg 22 


1297 


6.70 


0.78, - 


■1.37, 1.26 


-0.84, 


-4.39,0.43 


*CPD-46 8221 


2074573 


97 


1 


Hogg 22 


1297 


6.70 


-0.45,- 


-5.19, 1.32 


-0.84, 


-4.39, 0.43 


HD 150958 


2074576 


93 


84 


Hogg 22 


1297 


6.70 


-0.95,- 


-3.67,1.20 


-0.84, 


-4.39,0.43 


HD 153919 


1877055 





100 


NGC 6281 


494 


8.51 


1.42, 


4.84, 1.04 


-2.96, 


-3.75,0.23 


HD 156041 


1878768 


98 


63 


Bochum 13 


1077 


7.08 


0.03, - 


-0.74, 1.70 


-0.28, 


-1.20,0.30 


HD 156154 


1878846 


82 


100 


Bochum 13 


1077 


7.08 


-0.88,- 


-2.53, 1.64 


-0.28, 


-1.20,0.30 


*HD 164794 


1593528 


69 


100 


NGC 6530 


1322 


6.67 


1.92, - 


-0.40,1.09 


2.01, 


-1.81,0.51 



Notes. For each star, we give in the first two columns the object name and the corresponding ASCC catalogue number. The 
kinematic and photometric probabilities for cluster membership presented in Cols. 3 and 4 were calculated according to the 
procedures described by Kharchenko et al. 1)2004^ . The cluster names, the distances, and the ages are presented in the next three 
columns. The proper motions and their errors for stars and clusters are given in the last columns. All cluster data are taken from 
Kharchenko et al. (|2005ap . Stars with detected magnetic fields are marked with an asterisk in the first column. 



ity measurements) is assigned to several magnetic massive 
stars, such as (Oph (BlaauwT^), HD 93521 fGies fTW|) . 
HD 1 57857 (Schilbach & Roeser,2008j), HD 57682 (Comeron 
et al. nM5|) . and the thre e Of?p stars HD 108, HD 148937, 
HD 191612 (Hubrig et al. I2011dp . It is known that some of 
the massive stars located far from star clusters and star- 
forming regions are runavifays. They were likely formed in 
embedded clusters and then ejected into the field because 
of dynamical few-body interactions or binary-supernova ex- 
plosions (Eldridge et al. I201ip . Although the group of field 
O stars whose runaway status is difficult to prove does ex- 
ist, the recent study of Gvaramadze et al. (|2012p showed 
no significant evidence in support of the in-situ proposal 
for the origin of isolated massive stars. 

In some cases, even if the star is considered as an open 
cluster member, a careful study of their kinematic charac- 
teristics indicates that it is not fully co-moving with the 
other cluster members. As an example, the space motion 
of the magnetic 0-type star 0^ Ori C was studied by van 
Altena et al. (|1988p . who reported that 0"^ OriC is moving 
at 4.8±0.5kms~"'^ towards position angle 142° and that this 
velocity is significantly larger than the dispersion value of 
1.5±0.5kms~^ found for the other cluster members. The 
results of the radial velocity study of Stahl et al. (|2008p 
indicate that this star is moving rapidly away from the 
Orion Molecular Cloud and its host cluster. Another exam- 
ple is the kinematical status of the magnetic 0-type star 
HD 328856, which is on the basis of the photometric mem- 
bership probability a member of the compact open cluster 
Hogg 22 in the Ara region at an age of 5Myr and a dis- 
tance of about 1300 pc. However, its proper motions indi- 
cate that HD 328856 is not fully co-moving with the other 
cluster members, deviating from the cluster mean proper 
motion by ~2cr (Hubrig et al. 2011b). 

Among our studied runaway stars, HD 93521 and C, Oph 
are magnetic fast rotators with (3 Cep-type pulsations. 
These stars are thus particularly interesting to calibrate 
stellar rotation models including magnetic fields in the hot 
upper part of the Hertzsprung- Russell diagram. Variability 
linked to stellar pulsation has never been searched for or 
found in any of our other sample targets. In fact, there 



are only a few 0-type stars for which convincing evi- 
dence of non-radial pulsations with periods of a few hours 
has been detected from ground-based observations. Besides 
HD 93521 and COpli: there are ^ Per and ACep (de Jong 
et al. I1999p . For such stars, the photometric amplitudes 
are below mmag level, which are very difficult to detect 
from the ground. It may explain why only few pulsating 
0-type stars have been identified so far and mostly spec- 
troscopically. Recently, the very high-precision photome- 
try of the CoRoT satellite led to the detection of stellar 
pulsation in three other 0-type objects that belong to the 
young open cluster NGC 2244 inside the Rosette nebula and 
to the surrounding association MonOB2. /3 Cep-like pulsa- 
tion modes have been found in HD 46202 (Briquet et al. 
ECTTTl) and in HD 47129 (Mahy et al. lWTj) while unexpected 
modes of stochastic nature have been detected in HD 46149 
(Degroote et al. I2010P . Efforts to discover magnetic massive 
pulsators are worthwhile since their study provides a unique 
opportunity to probe the impact of a magnetic field on the 
physics of mixing inside stellar interiors. 



4.3. Indirect B-field indicators 

One of the major characteristics of magnetic Ap and Bp 
stars is the spectral variability due to the presence of chem- 
ical spots showing a certain symmetry with respect to the 
magnetic field geometry. Although spectral variability is 
present in all discovered magnetic massive stars, the inter- 
pretation of line profile variations seems not to be as simple 
as for Ap and Bp stars: the behaviour of line profiles of dif- 
ferent ions of a given element does not clearly correspond to 
the abundance variations. Synthetic spectra calculations of 
the line profile variations in massive stars are furthermore 
more troublesome since different spectral lines are to a dif- 
ferent degree sensitive to the stellar wind and to non-LTE 
effects. In any case, a detection of line profile variations 
can certainly be used as an indirect indicator for the possi- 
ble presence of a magnetic field. Another aspect that could 
hint at the presence of a magnetic field in massive stars 
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is the fact that a number of magnetic OB-type stars show 
nitrogen enrichment (e.g. Morel et ah I2008|) . 

It is not clear yet how to use the X-ray emission as 
an indirect indicator for the presence of magnetic fields. 
Anomalous X-ray emission can be expected from a mas- 
sive star with a magnetic field strong enough to regulate 
the stellar wind motion. The relative importance of mag- 
netic fields in gases is described by the plasma- /3p param- 
eter with /3p = 87rp/i?^, where p is the gas pressure. The 
gas is magnetically dominated when /3p < 1. For supersonic 
flows, such as stellar winds, the ram pressure exceeds the 
gas pressure and the dynamical importance of a magnetic 
field is defined by the ratio of wind kinetic to magnetic en- 
ergy density: when this ratio is less than one, the magnetic 
field dominates the wind bulk motion. 

Oskinova et al. (120111) investigated the X-ray emission 
and wind properties among the complete sample of known 
magnetic early B-type stars with available X-ray observa- 
tions. The authors have found that the X-ray emission from 
magnetic B-type stars is diverse. While some stars display 
hard variable X-rays, others are rather soft sources. The X- 
ray luminosities differ among otherwise similar Bp stars by 
more than two orders of magnitude. Among magnetic O- 
type stars, the situation is similar. While some of them are 
bright, hard, variable X-ray emission as predicted by the 
magnetically confined wind-shock scenario (e.g. Schulz et 
al. I2000p . others are relatively soft and faint X-ray sources 
(Naze et aL l^OTOI . 

Thus, the analysis of X-ray observations of magnetic OB 
stars shows that some of them are lacking strong and hard 
X-ray emission, or even are X-ray dark. Therefore, while a 
strong hard X-ray emission serves as a good indicator for 
the potential presence of a strong stellar magnetic field, X- 
ray strong stars represent only a sub-class of all magnetic 
massive stars. 

5. Outlook 

What is the origin of the magnetic fields in massive stars? It 
is quite possible that massive 0-type stars behave similar 
to slowly rotating, chemically peculiar, magnetic Ap and 
Bp stars, which amount for 10 to 20% of main-sequence A 
and B stars. It was suggested that their magnetic fields are 
inherited from their parent molecular clouds (Moss i2Q03l) . 
Alternatively, from analysing possible paths for the forma- 
tion of Ap and Bp stars based on modern theories for the 
evolution of single and binary stars, Tutukov & Fedorova 
(|2010p suggested that merging of close binaries is their main 
formation channel. 

Indeed, the evolution of a high fraction of massive stars 
is affected by environmental effects. Sana et al. (| 2012|) find 
two thirds of all Galactic O stars to be members of close 
binary systems. Mass transfer or stellar merger may rejuve- 
nate the mass gaining star while strong differential rotation 
may cause the creation of a magnetic field (Langer r2012|l . In 
dense star clusters, close multibody interactions may cause 
single and binary stars to be scattered out of the star cluster 
(e.g. Leonard & Duncan 1990). 

Clearly, these various ideas bear significantly different 
predictions for the appearance of magnetic main sequence 
stars inside and outside of star clusters. This means that 
these predictions are well testable, and we are confident 
that the question of the origin of the magnetic fields in 
massive main sequence stars will be answered in the near 



future. As described in Sect. SI first trends are emerging 
from the scarce data at hand. However, at the current stage, 
we clearly need more observations before firm conclusions 
can be drawn. 
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